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Abstract. We present experimental data on a one-dimensional superconducting 
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c/5 ' metamaterial that is tunable over a broad frequency band. The basic building block 

of this magnetic thin-film medium is a single-junction (rf-) superconducting quantum 

interference device (SQUID). Due to the nonlinear inductance of such an element, its 

resonance frequency is tunable in situ by applying a dc magnetic field. We demonstrate 

'"O . that this results in tunable effective parameters of our metamaterial consisting of 

^ ' 54 SQUIDs. In order to obtain the effective magnetic permeability /ir,cff from the 

fi ■ measured data, we employ a technique that uses only the complex transmission 

coefficient 6*21. 
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1. Introduction 

Losses and the strong limitation to a narrow frequency band are the main challenges 
when designing metamaterials that are made of conventional resonant structures. It 
has been shown that losses in metamaterials, working in and below the THz frequency 
range, can be greatly reduced when metallic structures are replaced by superconducting 
ones [I]. Additionally, superconducting meta-atoms exhibit an intrinsic tunability of 
their resonance frequency by magnetic field and temperature [21 [3]. However, in both 
cases, the tunability arises from a suppression of the superconducting order parameter, 
i.e. the density of Cooper pairs. Thus, by tuning the resonance frequency, the quality 
factor of the resonance is changed as well. 

In this work, we demostrate a one-dimensional metamaterial that employs 
superconducting quantum interference devices (SQUIDs) as meta-atoms, based on a 
theoretical idea introduced and further investigated in |31 15] and [6]. The SQUID can 
be considered as a split ring resonator (SRR) [7j that includes a Josephson junction. 
The tunability of its resonance frequency arises from the inductance of the Josephson 
junction that is tunable by a very weak dc magnetic field and does not come at the 
cost of suppression of superconductivity. The magnetic field tunability of the resonance 
frequency of such a SQUID meta-atom has been experimentally verified in [8]. This 
reference serves as basis for the microwave properties of a single SQUID. 

2. The rf-SQUID 

The basic building block of the metamaterial used in this work is a single junction 
(rf-)SQUID. Such an rf-SQUID consists of a superconducting loop interrupted by a 
Josephson junction (cf. figure[T](a)). The junction is indicated by a red cross. 



(a) 



-Dext 




Josephson 
junction 



shunt 



Figure 1. (a) Sketch of an rf-SQUID. The red cross symbolizes the Josephson junction, 
(b) Electric equivalent circuit for an rf-SQUID in the small signal approximation. £goo 
is the geometric inductance of the SQUID loop. The red circle indicates the electric 
circuit model for the junction. R represents the resistance due to a quasiparticle 
current, ij is the Josphehson junction and C stands for the capacitance between the 
superconducting electrodes, (c) Optical micrograph of the rf-SQUID. 



Within a small signal approximation, the Josephson junction can be considered as a 
nonlinear inductor [HI [ID]- When placed inside a superconducting loop, this so-called 
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Josephson inductance Lj is tunable by a magnetic field. In addition to Lj, the geometric 
inductance of the loop Lgeo contributes to the total inductance Ltot of the rf-SQUID. 
The full equivalent electric circuit is depicted in figureHJ^b). The red circle marks the 
electric circuit analogue of the Josephson junction for which the resistively capacitively 
shunted junction model is used [lO]. Like the SRR, the rf-SQUID can be interpreted as 
an LC-oscillator. Unlike the SRR however, the total inductance and thus the resonance 
frequency of the rf-SQUID is tunable, assuming the ac magnetic field component is 
small. 

Figure[l][c) shows an optical micrograph of the single SQUID. The SQUID and its 
junction are fabricated using a Nb/AlOx/Nb trilayer process. The Josephson junction 
is circular with a diameter of 1.6 yum, its critical current Ic = 1.8 fiA. From this value, 
the zero field Josephson inductance is calculated to be Lj = 183 pH. This value is 
approximately twice as large as the geometric inductance of the loop Lgeo = 82.5 pH. 
Thus, the rf-SQUID considered in this work is nonhysteretic [9], [10]. The junction is 
shunted with an additional parallel plate capacitor with a capacitance Cshunt = 2.0 pF 
which is two orders of magnitude larger than the intrinsic capacitance of the Josephson 
junction. Due to this shunt capacitor, the resonance frequency of the rf-SQUID is 
reduced and tunable between approximately 9 GHz and 15 GHz. 

3. The SQUID Metamaterial 

Two one-dimensional arrays of 27 rf-SQUIDs are placed inside the gaps of a coplanar 
waveguide (CPW) as shown in figure[2](a). The electromagnetic wave travels along the 
waveguide. In addition to the microwave signal, a dc current is applied along the central 
conductor, creating the dc magnetic field used for tuning the resonance frequency. Due 
to the waveguide geometry, the SQUIDs are oriented with their loop area perpendicular 
to the magnetic field. 

The pitch between neighboring SQUIDs in this sample is 92 fira, which is twice the 
width of the single SQUID and much larger than the distance between each SQUID 
and the central conductor of the waveguide. Therefore, the inductive coupling between 
adjacent SQUIDs is approximately one order of magnitude smaller than the coupling to 
the CPW and can be neglected. 

The CPW is connected to a vector network analyzer (VNA). The full experimental setup 
is depicted in figure[2](b). Rigorous magnetic shielding (not shown in the picture) proved 
to be crucial to protect the sample from stray magnetic fields. The sample inside the 
cryoperm magnetic shield, part of the attenuation and the amplifier are placed at liquid 
helium temperature, T = 4.2 K. 

4. Experimental Results 

We measure the complex transmission through the CPW (5*21) as a function of 
frequency z/ and magnetic flux $co- The microwave power at the sample is approx- 
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Figure 2. (a) Optical micrograph of part of tlie CPW containing a chain of rf-SQUIDs 
in each gap. (b) Measurement setup including the vector network analyzer (VNA), 
the biasTees, attentuation and cryogenic amplifier. 



iinately P ~ —90 dBin, including losses in the coaxial cables. The calibration of the 
measurement is done by applying a flux of $eo = ^cai = *^o/2. At this flux value the 
resonance frequency is shifted to its lowest possible value, which lies between 9 and 
10 GHz. The built-in "thru" calibration function of the VNA is used to substract the 
corresponding reference data from the rest of the measurement (see also supplimentary 
material). The resulting transmission magnitude for such a measurement is presented 
in figureEl For clarity, only data above 10 GHz are shown. 

Note that for the presentation of the data in this section, the more natural quantity 
of magnetic flux per SQUID $co is used instead of magnetic fleld. For our SQUIDs, a 
flux of one flux quantum $o = h/{2e) = 2.07 x 10^^^ Vs corresponds to a fleld of 1.2 /iT. 



16 



10 



1 1 1 r 



15 


- 




:;d 


N 






— 


ffi 14 
O 


- 


I 


-6 CN 

-12 


=^ 13 


_ 


>. 




o 




S 12 


- 







cr 




a 11 


- 


14-1 










/ 



>. 



_I L. 



-0.6 -0.4 -0.2 0.2 

external flux $eo/^o 



0.4 



0.6 



Figure 3. Measured transmission magnitude 15211 depending on frequency and 
magnetic flux $co/'&o- 



The resonance dip and its tunability with magnetic flux is clearly visible (low 
transmission is shown in red). It is the same behaviour that was already observed 
in [8] for a single rf-SQUID. Here, however, the resonance is much more pronounced 
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due to the increased number of SQUIDs. It should also be noted, that the shape and, 
more importantly, the depth of the resonance do not change considerably over a wide 
frequency band between 13 GHz and 14.5 GHz. The distortions in the onset and shape 
of the resonance dip, for example at 13.6 GHz and at 14.5 GHz, are due to parasitic 
sample holder resonances that couple to the SQUIDs. They are also faintly visible as 
flux independent stripes far away from resonance. 

The few separate lines, that are shifted in the flux axis with respect to the main 
resonance line, are lines from a few indivual SQUIDs. Their magnetic flux bias is 
different from the flux bias of the main part of the SQUID array. This is most 
probably due to trapped Abrikosov vortices either in the central conductor or in the 
superconducting material of the SQUIDs themselves. 
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Figure 4. (a) Transmission magnitude (top) and phase (bottom) at an external flux of 
$eO = — 0.185$o- (b) Real (top) and imaginary (bottom) part of the effective magnetic 
permeability /ir,eff calculated from the transmission data shown in (a). It should be 
noted that we employ /ir.off = Re(/Xi-^off) + iIm(/ir,off), opposite to the commonly used 
definition fj, = fi' — i fj," . 



FigurelD^a) shows a zoom to a smaller frequency range for a fixed fiiox value of external 
flux $eo = — 0.185$o- It allows a more detailed look at transmission magnitude and 
phase around the resonance frequency of z/q = 13.88 GHz. We observe that the resonance 
dip is indeed one collective resonance dip, no single SQUID lines are distiguishable. 
The steep rise of the phase dependence underlines the collective behaviour. However, 
as expected from the weak coupling between adjacent SQUIDs, there is no mutual 
synchronization. The quality factor of this resonance is QcoUective = 100. When we 
compare this value to the quality factor of a single SQUID line (we used one of the 
weak stray lines mentioned earlier), the quality factor of the single SQUID is more than 
twice as large Qsingie = 215. A detailed investigation of the quality factors, however, is 
hindered by the very shallow resonance of the single SQUID which is less than IdB. 
By using a method described in the next section, we are able retrieve the effective 
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magnetic permeability /Xr.efi from the measured transmission data depicted in figureO 
In figurelU^b), real and imaginary part of the effective magnetic permeability /ir,eff 
are plotted for the same frequency interval and flux value as the transmission data 
in FigJD^a). The magnetic permeability shows a frequency dependence typical for a 
metamaterial consisting of resonant elements [H]. Upon approaching the resonance 
frequency, the real part of /ir,eff increases from unity to almost 2. When reaching z/g 
it drops drastically to values below zero, which is due to the change of phase between 
incoming signal and magnetic response of the SQUIDs. From there it increases again 
to unity. The comparatively slow increase of Re(/ir_eff) between 12.5 GHz and 13.5 GHz 
as well as the small dip slightly above 13.5 GHz is due to a sample holder resonance at 
13.6 GHz which couples to the SQUIDs and affects the resonance. The deviation from 
zero of the imaginary part of /ir,cff reflects the increased losses at resonance. 
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Figure 5. (a) Frequency and flux dependent real part of the effective magnetic 
permeability Re(/ir,off) calculated from the transmission data shown in figure[31 The 
black dashed line indicates the cut shown in (b). (b) Flux dependence of Re(/ii.^off) at 
a frequency v = 13.83. Note that only negative flux values are shown. 



The full frequency and flux dependence of the real part of the effective magnetic 
permeability is shown in colour scale in figure[5]^a). The frequency dependent behaviour 
of Re(/ii.^eff) described above, extents qualitatively over the full frequency and flux range 
of the resonance curve displayed in figureO The data show clearly how the SQUID 
metamaterial enables us to reach any desired value of /ir,eff for a given frequency. The 
minimum and maximum achieved values are Re(/ii.^efT)min = —2 and Re(/ir,eflF)max = 
3, respectively. Unfortunately, due to sample holder resonances, especially around 
13.6 GHz and close to the maximum of the curve at approximately 14.5 GHz, these 
values cannot be reached for every measured frequency. Figure|5]^b) shows the real part 
of the magnetic permeability for a fixed frequency of z/ = 13.83 GHz, highlighted by the 
black dashed line in figure[5]^a). For this frequency, Re(/ir,eff) can be tuned to any value 
between -1.5 and +2. 
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5. Effective Parameter Retrieval 

To retrieve the effective, relative permeability from 5*21 for data presented in figuresHJ^b) 
and O we treat the system as a transmission line problem; the circuit diagram is shown 
in figure[6]^a). L' and C are the respective characteristic inductance and capacitance 
per unit length of transmission line. Then we simplify the unit cell by projecting the 
influence of the inductively coupled SQUIDs onto a change of the effective magnetic 
permeability seen by the transmission line and thus the line inductance (cf. Fig. Et^b)). 
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Figure 6. Electrical model of the circuit for one SQUID per unit cell. The unit 
cells in both pictures are indicated by a dashed box. (a) The CPW is modeled as a 
transmission line. The SQUIDs couple to the line via a mutual inductance only, (b) 
Further reduction of the system by including the SQUID influence into the relative, 
effective permeability ^r.off seen by the line. 



Using this model, it is a straight forward task to calculate the impedance matrix Z and 
scattering matrix S for the chain by cascading the unit cell's ABCD matrix. 
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Here, A = A {fir,cs) and A are the ABCD matrices of the unit cell and the chain of N unit 
cells, respectively. Zq is the port impedance on both sides. Since the transmission line 
properties L' and C of the CPW are known, the only free parameter is the effective, 
relative permeability /ir,eff- The reverse process, however, (calculating /ir^efr from the 
scattering matrix) is more complicated because it involves choosing the correct root 
when solving a system of coupled nonlinear equations (see supplementary material). 
Similar to other well known methods |12], this can be done by comparing the measured 
to the calculated phase delay. It should be noted that this retrieval method allows the 
reconstruction of /ir,eff using only one element of the scattering matrix which simplifies 
the experimental process significantly. 



A one- dimensional tunable magnetic metamaterial 8 

6. Conclusion 

In this work, we liave presented experimental results on a superconducting thin-film 
metamaterial that is tunable in situ by a weak dc magnetic field. As basic building 
blocks, we employed rf-SQUIDs embedded inside the gap of a coplanar waveguide. We 
have demonstrated a tunable, collective resonance of almost all 54 SQUIDs. Changing 
the resonance frequency hardly affects the quality of the resonance in wide frequency 
band. However, sample holder resonances distort the picture for some frequencies. 
Using only the transmission data 5*21, we were able to extract the effective magnetic 
permeability from the measured data. These results show that we have created a 
material with tunable effective magnetic permeability in a frequency range between 
10 GHz and 14.5 GHz. Further experiments will include an improved sample layout and 
environment in order to obtain results with higher quality factors and less parasitic 
resonances. 
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en ■ Supplementary Material 

o. 

^ ' 1. Retrieval of yUreff 

5h ' 

^ ! In order to retrieve tlie effective parameter /ir,eff mentioned in Sections, we describe the 

f>^ ! loaded waveguide as a lossless transmission line (see Fig. 6). The ABCD matrix of one 

<N ; unit cell [T] is 

1: -4=1 "°T "■: I- (1) 

^ ■ Here, Zl = iuj^r,esL'l is the impedance of the inductance and Zc = l/{iujC'l) is the 

^ ■ impedance of the capacitance for a unit cell of length /. L' and C are the inductance 



and capacitance per unit length of transmission line, respectively. The eigenvalues 



'l - TT^ ' '2 - ^-^ (2) 

Ch ■ and eigenvectors 

of the ABCD matrix can be used to rewrite it in the form 
^: A = C-{^' ^\-C-' (4) 

en . V /2 ; ^ ' 

Tt ! with C = [ 61,62 ) . (5) 

O 

cn ' This simplifies the total ABCD matrix A when cascading N unit cells: 



/f 
/^ 



A = A-=C.\ \ ,; ].C-\ (6) 



X 

c^ . Using this equation and the well-known relation [T] of A to the scattering matrix S^ we 

construct a system of four coupled, nonlinear equations 

{h-i)q -{h-i)i? {i?-i^)zi. 

h-h h-h I _ (7) 

(ii~i)i^h-{ii-i)i^-ah-i)h-h+i)i^ {h-i)i^-ifi2+i^ I ^^ 

{h-l2)Z-L h-h 

(Sll+l)S22-Sl2S21-Sli-l {(Sll+l)S22-Sl2S21+Sll + l)Zo 

2 S21 2 521 

(Sll — l)S22—Sl2S21—Sll+l (5ii— 1)522 "■5'l2 5*21+511— 1 

2 521^0 2 521 

for the four variables 6*11, 5*12, 5*22, and /ij.,eff since 5*21 is the measured quantity and 
all other parameters (like L' and C) are known from either design considerations or 
simulations. 
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In general, this system of equations is not single- valued. There are several ways to sort 
out implausible solutions like requiring reciprocity, energy conservation or comparing 
the measured to the calculated phase delay. 

2. Data calibration & reference planes 

In order to retrieve the effective, relative permeability from 5*21 using the method 
described above, special care has to be taken to choose the correct calibration technique. 
In this approach, the reference planes of the calibration have to be just before and after 
the array of unit cells. In a real experiment, the best set of reference planes achievable 
for a full calibration of such a measurement is located at the microwave connectors 
closest to the sample at cryogenic temperatures. This, however, does not only require a 
more complex experimental setup [2], it is also insufficient for our retrieval method. 

To solve this issue, we combine measurement and simulation to prepare the data for 
the retrieval algorithm. First of all, we divide the experimental setup into three parts, 
each of which can be described individually by a two-port scattering matrix: 

5*™ describes the input part of the setup, from port 1 of the network analyzer 

to the beginning of the SQUID loaded section of the CPW. 
S'^^ describes the SQUID loaded section of the CPW. 

^out describes the output part of the setup from the end of the SQUID loaded 
section of the CPW to port 2 of the network analyzer. 

We measure the total transmission from VNA port 1 to port 2 which can be writteiiji 

as 

Cin cstl cout 
otot ^21'-'21^21 /Q\ 

f Cin Cstl Cstl f Cin Cstl i A CstA Cout _|_ Cin Cstl i ' ^ ' 

\^^22^12^21 \^^22^11 ^J'J22J'-'ll "r '-'22^11 -■- 

If we assume that the waveguides leading to the sample are well matched (i.e. no 
reflections at the waveguide connections to the sample: 5*22 = 3°^^ = 0), the expression 
simplifies to 

Ctot cin cstl cout /q\ 

•^21 — '^21^21^21 • WJ 

Clearly, all of the above matrix elements are frequency dependent. 5*1*^ however, also 
depends on the dc magnetic field. It can be decomposed into two factors by assuming 

Sfi (w) is the transmission through the loaded waveguide section without the 
effect of the SQUID resonance. Consequently, this factor does not 
depend on magnetic field. 
a {u, $00 ) describes the change in the transmission through the loaded waveguide 
due to the SQUID resonance. Therefore, this factor is frequency and 
field dependent. 

I By multiplying the transfer matrices corresponding to these scattering matrices. [1] 
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The total transmission then reads 

Slf {iO, $eo) = 4l (^) Sfl (^) « (^, '^eo) 5°^ (^) • (10) 

The goal of the calibration is to extract the term 5"!^' from the measured Sg"* since 
the former is the quantity needed for the effective /ir,eff retrieval. If we restrict our 
investigation to a limited frequency range bounded by Ujam and Wmax, "we can usually 
find a value of the flux $oo = "^cai for which a {u, $cai) ~ 1- Thus, by dividing all 
the measured data in the specified frequency range by the corresponding value at the 
calibration fluji§|, we can extract a: 

otot.cal / ^ N '-'21 (^; 'I'co) ^ / ^ \ /li\ 

S^i' {U, 'I'eo) = ,,, . ^a{u, 'I'eo) (11) 

i>21 {^, "Pcalj 

As a last step, we have to reconstruct Sfi from the calibrated data S21 ''^^ . Therefore, we 
simulate the transmission through the loaded waveguide without the SQUID resonance 
Sfi (u) and multiply it by the calibrated data: 

ristl, reconstructed / /r \ ritot,cal / ,t, \ ristl, simulated / \ /i o\ 

S21 {(^,^cO) = S21' {uj,^co)- S21' [UJ) (12) 

This result can then be used in the algorithm outlined in the last section. It should be 
stressed, however, that this method is only valid if the effect of the SQUID resonance 
at the calibration flux is negligible in the frequency range of interest. 
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§ Choosing this "cahbration" flux effectively limits the validity of the measured data to the mentioned 
frequency range. 



